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1. Introduction 

The reconstruction of tiie properties of an object from data obtained in a scattering 
experiment is a problem that arises in all branches of physics. Here we consider the 
scattering of light from an object placed inside a diffuse medium. This poses problems 
because of the many further scattering events between object and observer; in everyday 
language, "the object cannot be seen" . However, with the help of sufficiently intense light 
sources, and careful and sensitive experimental techniques the small perturbations on the 
diffuse light due to the embedded object become detectable. For example experiments 
have been performed by den Outer, Nieuwenhuizen and Lagendijk 0. The technological 
applications to non-invasive methods of detecting tumors in human flesh are active areas 
of research. 

The significant feature of this work is that we go beyond the diffusion approximation 
discussed in and make analytic calculations based on mesoscopic theory. This means 
that although the scattering length is much longer than the wavelength of the light {kl ^ 
1), we do not require that the relative size of the objects be much larger than the scattering 
length. Furthermore, we are able to impose boundary conditions based on the microscopic 
theory, which cannot be done in the diffusion approximation. We start with the simplest 
case, namely a point object. Our main result is for a small spherical object and we expect 
this result to be useful in providing an independent check of the numerical and other 
approximate approaches sometimes used in this field 0. 

We make several simplifying assumptions: most importantly we only consider scalar 
waves and always assume isotropic scattering. Our results are presented in the form of 
the effective charge and dipole moment of the multipole expansion of the solutions of the 
diffusion equation in the far field region. This allows the results to be adapted to whatever 
detection arrangement is desired. 

The paper starts with a description of the salient points of the mesoscopic multi- 
ple scattering formalism forming the basis of our analysis. The following sections, 3,4,5 
consider a series of progressively more realistic objects: a point scatterer, a collection of 
such scatterers and a small sphere subject to a variety of boundary conditions. A short 
conclusion summarises our results. 
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2. Mesoscopic Scattering Formalism 

We base our approach on the Schwarzchild Milne (SM) equation of radiation transfer 
theory. This equation may be obtained by integrating the full radiation transfer equation 
but here we prefer to regard it as a Bethe Saltpeter-like integral equation for the diffuse 
intensity J(r). We write it in the general form, 

J(r) = S{r) + J d^r'M{r,r')J{r'). (2.1) 

Where the source 5" and the kernel A4 take different forms depending on the geometry of 
the problem. In the ladder approximation the kernel is given by the square of the amplitude 
green functions: Ai = ^[^(r, r')p. A uniform density of scatterers leads to an amplitude 
green function Go{r) = e^^^^'^^'^^^'^ /Aur, where the scattering length / is related to the t- 
matrices of the individual scatterers and their density (n) by / = 4n/ntt. It is then simple 
to see that M.o = e~'^^^ / Anlr'^ (the SM equation is often stated for a slab geometry, and 
the relevant kernel is then obtained from TWq by integrating over perpendicular momenta). 
The Fourier transform of this kernel is M.Q{q) = taii~^{ql) / ql (= Ai[ql) in the notation 
described in appendix C), so the green function Qq for intensity propagation in the SM 
equation is Qoi^q) = 1/{1 — Ai{ql)). A fuller discussion of notation is provided in appendices 
A and B, and a general review of the formalism is given in 0]. 

In order to find the effect of an object embedded in the multiple scattering medium 
we require an expression for the appropriate SM kernel, which we base on a perturbative 
expansion of the amplitude green function. Firstly however, we consider exactly what 
quantity should be calculated to compare with measurements. 

2.1. Measurement 

We assume that system boundaries are distant from the localised region of the object. 
In this far field region, the diffuse intensity J(r) obeys the diffusion equation which is 
much simpler than the full SM equation. Since we are dealing with static quantities, the 
diffusion equation reduces to the Poisson equation. 

J(r) = 0, (2.2) 

In the absence of the object, the solution of equation (pT2|), subject to appropriate boundary 
conditions, is Jq{t). The total intensity in the presence of the object then reads 

J(r) = Jo(r)+(5J(r). (2.3) 
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The goal of the present work is to characterise the disturbance 5J{r) of the diffuse intensity, 
representing the diffuse image of the immersed object. We make use of the fact that the 
Poisson equation ( |2.2| ), allows electrostatic analogies, and represent the disturbance dJ{r) 
in terms of a multipole series. It turns out that only the first two of them, namely the 
charge q and the dipole p, play a role, while higher-order multipoles are negligible. Assume 
for a while that the turbid medium is infinite. The disturbance of the intensity far enough 
from an object placed at Tq assumes the form 



= - (2.4) 

|r — fqI |r — fqI'^ 



This expression, which is a solution of the diffusion equation ( pT2| ) , is also a solution of the 
full SM equation provided the distance from the object is large with respect to the mean 
free path 

For a symmetric object, the linearity and isotropy of the problem imply that the 
charge and dipole have the form 

Q = -QJo(ro), p = -PVJo(ro). (2.5) 

We follow the notation proposed by J.M. Luck, and call Q the capacitance of the object, 
and P the polarizability [||. These two numbers are intrinsic to the embedded object. 

The capacitance Q is non-zero (and in fact positive) only if the sphere absorbs light. 
For a non-absorbing object, the leading contribution is that of the induced dipole p. The 
corresponding polarizability P can be either positive or negative, depending on the nature 
of the object. Notice that Q has the dimension of a length, while P has the dimension of 
a volume. 

For real measurements we should bear in mind that the diffusion approximation will 
break down in the vicinity of the measurement apparatus. For the particular case of a thick 
slab geometry, which is often encountered experimentally ^ , we now recall the derivation 
of the diffuse image of a small object to show how the capacitance and polarizability can 
be measured. 

The sample is a thick slab (0 < x < L), with L ^ /, assumed to be infinite in the 
two transverse directions. We set p = (y^ + z^)^/^. The left side of the sample (x = 0) 
is subject to an incident plane wave, so that the solution of eq. (|2.2|) when the immersed 
object is absent reads 

Jo(r) = (2.6) 



The intensity T(y, z) transmitted through the sample, and emitted on the right side (x = L) 
at the point {y,z), is proportional to the derivative dJ{x = L,y, z)/dx. When there is 
no immersed object, eq. ( |2.6|) yields a uniform transmission Tq = Kl/L. The prefactor 
K describes the angular dependence of the transmitted intensity. It also depends on 
microscopic characteristics of the medium, such as any anisotropy of the scatterers. 

Let the immersed object have capacitance Q and polarizability P, and be located at 
X = xq, y = z = Q. Its charge and dipole then read 

? = -Qj(l-f), P=^£. (2-7) 

with X being the unit vector of the positive x-axis. 

Along the lines of ref. [Q, the total intensity can be determined by summing the 



free-space expression ( p.4| ) over an infinite double array of images (—00 < n < +00). We 
thus obtain 

x\ ^ / 1 1 



J(x) = j(l-|)+,5^ 




[(x - xo + 2nLY + p2j V2 ^1^^ + xo + 2nL)2 + p2j V2 
X — xo + 2nL ^ X + Xo + 2nL 



[(x - Xo + 2nL)2 + p2] 3/2 ^^^^ 2nL)2 + p2] 3/2 



The transmitted intensity T[y^ z) = T{p) then reads 

L — Xo + 2nL 




n [(L-xo + 2nL)2+p2]3/2 
2(L-xo + 2nL)2-p2 



(2.9) 



[{L -xo + 2nL)2 + p2] 
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Such characteristic profiles have been observed experimentally for the case of pencils (black 
cylinders) and glass fibers (transparant cylinders with index mismatch) by den Outer et 
al. 

In order to calculate the charge and dipole moment in terms of microscopic properties 



of the object we relate (|2.4|) to the far field behaviour of the full SM equation with the same 



incoming intensity. At first order, an object placed in the medium modifies the free kernel 



of the SM equation by adding a vertex term Ai = Aio + V, which for smaU momenta we 
can writeB as, 

V{p,p') = ^{-Q + Pp.p'). (2.10) 

The parameters Q and P are precisely the ones that appear in the multipole expansion 
( p.4|) . To see this, insert the expansion of the kernel, Ai = A4o + V, into the the SM 
equation and solve perturbatively for the far field intensity J(r). This can be done using 
the diffusion limit of the free intensity green function, Qo{r) oc 1/r. We emphasise the 
pertubative nature of this general result, and indeed it fails in the case of a large spherical 
object. This will be the reason for later restricting our attention to a small sphere. 

In the rest of the paper we will define the vertex more precisely and develop a dia- 
grammatic perturbation expansion of the SM equation to calculate it both at leading and 
higher order. Only the small momentum behavior of the vertex is important for our goal 
of calculating Q and P, and we shall find that in this case the higher order corrections can 
sometimes be resummed. 

3. Additional Point-Like Scatterer 

The effect of an object consisting of a single point-like scatterer is a worthwhile exercise 
because it allows us to demonstrate the method in a simple setting. Although the higher 
order terms can be formally calculated, they are not reliable since we encounter ill-defined 
expressions arising from the point-like limit which take us out of the realm of validity of 
the mesoscopic theory. 

3.1. Leading Order Calculation 

Starting from the amplitude green function for a single extra scatterer with t-matrix 
te, at the origin, 

^(r, r') = Godr - r'|) + Go{r%Goir'). (3.1) 
The full SM kernel is given by, 

4-7r Att 

M{r, r') =Mo + 5M = —\G\^ = —\Goi\r - r'|) + Go{r%Go{r')\^ (3.2) 



^ We define the Fourier transform so that the first momentum is incoming and the second 
outgoing. 
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Fig. 1: Diagrams contributing to the bare vertex Vi. 



The diagrammatic interpretation is obvious and the vertex ( p.lO| ) contains the terms 
shown in figure 1. AU terms wiU be the same order due to the optical theorem. If we 
evaluate the diagrams in momentum space we may use the formulae in appendix B to 
obtain the mesoscopic limit, 

Vi(p, p') = dM{P. P ) = I (^-^AMA^ip'l) + ^^^^^^2(P/, P'/)^ . (3.3) 

Now we use the optical theorem as applied to the additional scatterer. We allow 
this extra scatterer to have a complex refractive index and introduce the albedo which 
appears in the optical theorem: 

(yei = -r—=ae — ; — (3.4) 
47r k 

Where a^i is the elastic contribution to the cross section for the additional scatterer. 
The final expression for the bare vertex is, 

yi(p,p') = la,i(^A^{pl)A^{p'l) - 1a2(p/,p'0) 



(3.5) 



By using the small momentum expansions of the Ai functions given in appendix C, we are 
able to identify the charge and dipole moment defined by equation ( p.lOj) , 

^ CTg; 3 3 aabs 



I 4:71 A-K I /Q a\ 

I I (3-6) 

Pi = (Jeil — il/ae) = —(Jtotl 

471 47r 



For albedo Oe < 1, there is absorption and the charge is positive. 



3.2. Second Order Correction 

This expression for the vertex, Vi = 6Ai, is only the first term in the perturbative 
expansion of the SM equation. Higher order terms contain these vertices connected by 
intensity propagators in the form Vn = ViQqVi . . . ViQqVi. In general these vertices are 
not easy to evaluate because the intermediate loop integrals do not give rise to closed 
expressions. However, for the multipole expansion we only require the behaviour for small 
external momentum and can therefore approximate the first and last Vi(p, p'). In this 
case, the angular integration over the intermediate momenta simplifies the expressions 
remarkably so that the effect propagates though the series making it valid for each loop as 
shown below. 

For small incoming momentum, p, we can expand the bare vertex using the expansion 
of A2 in appendix C. This leaves Vi(p, p') written in terms of Ai{p') and explicit factors 
of momenta: p^, p'^ and p.p'. Let us work to first order in p in which case, 

^i(p,p') -^QiA,{p'l) +p.p'47rl'P^As{p'l) (3.7) 

where Qi and Pi are given above in and we have employed the useful combination 
As{q) = il-Ai{q))/q^. 
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Fig. 2: The second order vertex V2, is composed of two first order vertices. 



Now consider the second order contribution shown in figure 2, 

^2(P,P') = I -0^ ^i(p,q)^o(?)^i(q,p') (3.8) 

where Qq is the free propagator as dicussed in section 2. Using the approximate expression 
( pi.7| ) in the first factor of Vi the integrals become tractable. Performing the angular part 
of the integral using the formulae given in appendix C we find that the general form of V 
is preserved, 

V2{p, P') = -^Q2Ai{p'l) + p.p'4nfP2As{p'l) (3.9) 
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where the exphcit expressions for Q2 and P2 are, 



(^ei f ^ ,\ ^ 47r ^ ^2 



Q2 = ^ - - 1 SiQi = -^B,Qi 

^ (3.10) 

P2 = -(Tel ( ] B2P1 = -AttB^P! 



a, 



The remaining integrations over the magnitude of the momenta are contained in the di- 
mensionless factors Bi and B2, where in the second expressions we have used the exphcit 
form of Goiq). 



r r 1 r 

Bi = ^ J dq q^Aj{ql)goiq) = 



u ^'-^'^'^^ (3-11) 
^^^2^ J "^"^ ^^'^siQl)Goiq) = ^ dq (1 - A^{q)) 

Here we encounter a problem due to the point-hke hmit. The integrals diverge linearly, so 
we have introduced a high momentum regulator A which we expect to be of the order of 
the inverse size of the scatterer, A ~ 1/R. For the point-like analysis above to be valid, 
the size of the extra scatterer must be much less than the wavelength, R << X. This 
means that very high momenta, up to the inverse scale of the scatterer, are encountered 
in the integral and in this range the mesoscopic approximation fails. The expressions for 
the second order terms are therefore formal. A similar divergence arises for the t-matrix of 
a point scatterer , where it was parametrised as a linear divergence depending on some 
property of the scatterer. The problem of divergencies will not appear in the following 
sections where we will consider extended objects. 

3.3. Higher Order Corrections 

Bearing in mind the caveats above, we proceed formally to illustrate how the higher 
order effects may be calculated. Since the form of the vertex remains the same as additional 
bubbles are added, it is clear that the full series can be summed to obtain an expression 
for the full V. The explicit expressions for Q and P are then, 

^_ Qi _ Qi 

P= ^ = ^ 

l + 3aJ{l/a,)B2 1 + 47rPiS2 

It is interesting to note that if the capacitance Q is negative (corresponding to amplification 
rather than absorption) it could be substantially enhanced by these higher order effects. 



8 



3.4- Maximally Crossed Diagrams 

Because we are dealing with the pecuhar situation of a point-hke object another issue 
becomes important, which although academic is interesting to discuss. Crossed diagrams 
generally play a sub-leading role in mesoscopic calculations, for example in the ladder 
propagator it can be seen that crossed diagrams contribute at higher order in l/kl. This 
conclusion must be reexamined for the point object and diagrams in real space show that 
the crossed diagrams can contribute at the same order as uncrossed ones. This can be 
understood as being due to the phase cancelling in the same way as it does for enhanced 
backscattering. These contributions only occur when maximally crossed series are inserted 
between two additional scatterers on the same line (there can be an extra scatter on the 
other line or not). In this case, each diagram contributes exactly as for the ladder sum. 
The previous calculation must be generalised so that crossed diagrams are not included 
between additional scatterers on different lines. This combinatoric problem may be solved 
using a matrix formalism described in Appendix D and the final result in the mesoscopic 
limit is that the charge and dipole in that case are as above but B gains a factor of 2. 

These calculations illustrate how the full sum can be calculated from an infinite series 
of diagrams. Unfortunately the results are flawed because the divergences due to the point 
nature of the object take the calculation out of the realm of validity of the mesoscopic 
approximation. In the next sections we will consider better defined systems. 

4. Many Extra Scatterers 

We now consider an object consisting of a collection of extra point-like scatterers 
localised in some region of size R. We first develop perturbation theory for the amplitude 
green functions which we then use in the kernel of the SM equation. 

4.1. Perturbative expansion 

Consider extra scatterers placed in a localised region with density profile ne(r). The 
perturbative expansion of the amplitude green function in terms of the individual extra 
scatterer t-matrix tg, resembles the expression (PTT|), 






+ . . . 




9 



Where, as usual, Gq is the green function in the presence of a constant density of back- 
ground scatterers. In momentum space we have, 

G(p, p') = S{p - p')G'o(p) + ne{p- p')teGo(p)Go(p') 

r fPo (4-2) 
+ J ^n,(p-q)ne(q-p')ieGo(p)G'o(q)G'o(p') + ... 

We use the above expression to find a form for the SM kernel. Here we must be careful to 
include an extra term in the definition of the kernel because the amplitude lines may end 
on an extra rather than background scatterer; to the factor 4n/l interpreted as n|tp, we 
must add the contribution of the extra scatterers ne(r)|ieP- 

2 \ |/-y/ /m2 



A^(r,r') = (^ + ne(r)|te|2^|G(r,r 



-M(p,p') = y/ (2l)3(fj)3 ^(q+p-q^+p^)^(q-q^) (^.s) 

In this problem, the vertex that relates the mesoscopic SM calculation to the multipole 
expansion ( p.4|) is not simply given by ^ = 6Ai as it was for the single point object. In the 
far-field region where the identification is made, there are no extra scatterers, so we must 
insist that the vertex only contain diagrams which end on a background rather than extra 
scatterer. This prescription is most easily viewed diagrammatically rather than in terms 
of the formal perturbation expansion of the SM equation and its effect is also to remove 
the asymmetry apparent in the kernel itself. 

Another difference with the point scatterer case is that the perturbation expansion is 
organised in powers of the extra density factor 7ie(r), and we are unable to relate different 
orders with the Optical theorem. The additional term included in (|4.3| ) saves the situation. 
Inserting the expansion of the green function and identifying the first order term A4 1 , the 
additional term also contributes as the last term below. 

-Mi(p,p) 



/ J (27r)3(27r)3 

tene(?-?')G'o(q)G'o(q')^o(q + p)5(q + p-q' -p') + c.c. (4.4) 



/d q 
^^^G'o(q)G'o(q + p) 
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This expression can be simplified using the functions Ai defined in appendix B and C. In 
order to obtain Vi from A^i, we make a small modification to the second term in order 
that it end on a background scatterer. 

^l(P,P') = (^^iteheip - p')-^A2{pl,p'l') + C.C^ + \te\'^ Ueip - p') ^Ai{pl) Ai{p' l') 

= n,{p - p')la,i (^A,{pl)A,{p'l') - ^^2(p^ p'l')^ 

(4.5) 

The second line follows from the optical theorem applied to the individual scatterers, as in 
the point scatterer case. We can immediately identify the charge and dipole at this order 
by taking the momenta small and using the definition ( [2.101 ). 



Pi = a J -^iVe l/tte - 6^(l/ae - 1) = N.Pr''' - 2Ar,p2QPOznt object 



47r 

Where A^e = ^e(O) is the total number of extra scatterers. The second term in Pi comes 
from expanding the density form factor for small momentum as he{p) ~ Ne{l — p'^ R'^) , and 
P is a measure of the size of the object, R^ = J r'^ne{r)dr. A more precise expression 
can be obtained if we assume a particular density profile, for example a Gaussian form. 
We have assumed a spherically symmetric density perturbation so nf,{p) is real. The 
similarity with the single scatterer result (|3.6|) is clear; the charge merely gets multiplied 
by the number of scatterers present whereas the dipole moment also picks up a term that 
depends on the size of the object. 

4-2. Second Order 

Two types of term occur at second order, those from second order in the expansion of 
the amplitude green function, and those from first order in the expansion of the amplitude 
green function but second order in the expansion of the SM equation. 

The first type of term can be neglected because for example one contribution to the 
charge is, 

^""itel' ! 7#!:,7^|ne(q-q')nGo(q)nGo(q')P (4.7) 



/ ' 7 (27r)3 (27r)3 

which corresponds to processes in which the amplitudes scatter off different extra scatterers 
so the phase cancelation kills them. This can be seen explicitly for a Gaussian profile, where 
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the contribution will be of order a^/l"^ x N'^/k'^E?. By comparing with the main second 
order contributions below we obtain a condition on the size of object for this analysis to 
remain valid, k'^B?{l - l/ae)^!^^"* > 1. 

We therefore only consider the term with the form V2 = ViQVi. Note that although 
this is clear diagrammatically, at the level of the formal expansion the correction we made 
to obtain Vi from M.i appears in ViQVi in just the right way to count the relevant terms 
in M.2-: which need not be evaluated separately. We attempt to use the same trick as for 
the point scatterer and therefore expand the bare vertex for small incoming momentum p, 

^i(p, pO ne(p - pO (^-^gr^* "''^^'A^ip'l) + p.p'47r/2pf°-* ''''''''' A^{p'l)^ 

(4.8) 

Now consider the first diagram, 

^2(P,P') = /(03 ViMOmiM) (4.9) 

Because of the presence of the density form factor, the angular integrals that kept the form 
of V the same at higher order and therefore allowed us to perform the full sum in the point 
scatterer case cannot be done. Let us instead do the second order calculation by setting 
both the external momenta small. The final expressions for Q2 and P2 are then, 

Q2 = ^iVe f 1 - —) B^Qi 
P2 = aelNe—B2Pi - 2R^Qi 

Where the factors B are now well-defined and do not need regulating because the den- 
sity acts as a smooth cutoff at momenta of order 1/R. In this regime the mesoscopic 
approximation is still valid so the approach is consistent. 



B 



1 f°° A'^(q) 

^'=2^^^i dqnl{q/l){l-A,{q)) 

This resolves the difficulties encountered in the point object case, but the discussion 
above shows that we cannot write simple expressions for terms higher than second order. 
Crossed diagrams are also suppressed in this case. 
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5. Spherical Object 

We consider tlie simplest example of an extended object immersed in the multiple 
scattering background: a sphere of radius R {kR ^ 1) subject to a variety of boundary 
conditions. The analogous problem in a free background was much studied in the early part 
of this century , and some of the results relevant to the present problem are collected in 
appendix E. 

The general result relating the effective charge to the fourier transformed vertex is 
pertubative about the free theory, but depending on the boundary conditions, the interior 
region of a finite size object may not be susceptible to such treatment. For example, for 
reflecting boundary conditions the full green function in the interior vanishes, which cannot 
be regarded as a small correction. A different approximation scheme, or an object with 
similar refractive index to the medium could be considered, but here we shall proceed by 
assuming that the sphere is small: R <^ I. In this case it can be argued that the correction 
from the interior of the sphere is smaller than the perturbative term we shall calculate. 
This limit is the opposite of the limit in which diffusion theory is valid. An alternative 
approach to this problem is to use radiative transfer theory and boundary conditions for 
the intensity field 0. 

5.1. Leading Order 

Employing spherical polar coordinates we write the amplitude Green function in terms 
of spherical harmonics and spherical Bessel functions, 

G(r, r') = J2 9i{r, r')YUO, <P)YUe\ <P') 

Im (5.1) 

gi{r, r') = g\'^\r, r) + aig\^\r, r') = ik[ji{kr<) + aih'l^\kr^)) h\^\kry) 

Where the notation r>, r< indicates the larger or smaller of the two radii r, r' . As usual, 
k has an imaginary part to take account of multiple scattering, but on the surface of 
the sphere where the boundary conditions determining the coefficients ai are imposed, 
the value of k must be reconsidered. In a fully self consistent approach, the imaginary 
part of k is fixed by the return green function G(r, r) 0], and can vary in space. This 
variation is only significant close to the surface of the extended object and in the mesoscopic 
approximation it can be ignored^, thus justifying the form ([5.11 ). At the surface however, 

^ In a one- dimensional version of the problem this can be an important effect. 
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where the boundary conditions are imposed, the self consistent form of k must be used. 
For example Dirichlet [D) and Neumann (A^) conditions which correspond to reflections, 
have k real at the boundary and lead to af^'^ = —(1 =F (— l)'e~^*'^^)/2. More general 
boundary conditions, mimicking the matching of internal and external solutions, can be 
imposed by introducing a surface impedance. 

In the same way as for the point-like object of section 3, the bare vertex of the SM 
equation is given by, 

^i(r,r') = y (|G(r,r')|^-|G('^)(r,r')p) (5.2) 

Where G*-'^^ is the free green function, given by the same expansion as ( |5.1|) with ai = 0. 
According to ( |2.10| ), the charge Q is related to the Fourier transform of this vertex at 
vanishing momenta. This pertubative result is clearly inappropriate in the interior of the 
sphere where G(r, r') vanishes, nonetheless corrections from the interior should not exceed 
0{R/l)^, and can be neglected in comparison with the leading term for small spheres. 
Using orthogonality properties of the spherical harmonics, the charge can be written as 
a sum of radial integrals over products of gf^\r,r') and g^^\r,r'). These integrals may 
be simplified using the asymptotic form of the Bessel functions valid since kR ^ 1. This 
asymptotic approximation is only good for angular momentum numbers / < Imax, where 
Imax = kR, which in any case forms a effective upper limit in the angular momentum sum 
since the a^'s become small for higher / (this approximation is standard in dealing with 
scattering from spheres). Dropping subleading terms in 1/kl we find. 

The combination of a^'s is recognisable from the classical amplitude scattering problem 
where a similar sum appears in the absorption cross section. Using this result, and only 
keeping the leading term in R/l, we find. 

Where aabs is the absorption cross section. It is reassuring that this coincides with the 
point-like result (|3.6|). In particular, for a black (totally absorbing) sphere one has aabs = 
ttR"^. This implies Qi = 3i?^/4Z, a result derived already in 0. 
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For the D and boundary conditions the charge must vanish since they correspond 
to reflection. There is no absorption and indeed each term in the sum ( p.3| ) vanishes when 
the exphcit forms of a^'^ are used. The calculation of aabs in the case of more general 
boundary conditions (for example a surface impedance) is a standard problem, the only 
subtlety arising from the value of k at the boundary. 

The dipole moment can be calculated in a similar way. We find the general result, 

(5.5) 

The first term in the sum is recognised as the total cross section otot = 27ri?^, but the 
second term is a combination that never appears in the amplitude problem. For boundary 
conditions based on a surface impedance we can show that this second term is proportional 
to the absorbtion cross section (and therefore vanishes for D or N boundary conditions). 
Taking only the leading order terms we may write the dipole as, 

= 7^ - ^ ^'-'^ 

The first term recovers the point object result in the limit R <^ I, whereas the second term 
is reminiscent of the similar term in the formula (^]^) for a collection of point scatterers. 
In case of a small refiecting but non-absorbing sphere we have atot = Stt-R^, and we find 
Pi = R^l/2. 

5.2. Higher order 

We may sum all the higher order corrections to the charge and dipole moment following 
the same method as for the point scatterer. In this case we do not have an explicit form 
for the Fourier transform of the vertex so intermediate formulae are unwieldy. We start 
with the form of the vertex for small incoming momentum in analogy to formula (|3.7|) . 

^i(p,p') -^QiAf\p'l,R/l)+p.p'A7,fP^A\^\p'l,R/l) (5.7) 

Where we have introduced functions A^^\p'l, R/l) and A^^\p'l, R/l), defined in appendix 
F. These functions become the point scatterer functions Ai, A^ in the limits of small R, 
but are better behaved at large momenta thus avoiding the problems encountered with the 
point object. 
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A calculation of the second order vertex at small incoming momentum shows that it 
retains the same form as ( |5.7| ). This calculation requires various integrals that are collected 
in appendix F. Proceeding by iteration the full sum can be performed exactly in the same 
way as for the point object. Our final result takes the same form as ( |3.12| ): 

n = 9l 



P 



3"^"^ (5.8) 
Pi ^ ' 



1 + AtiPiB2 



Where the integrals over the magnitude of the momentum are now given in terms of the 
modified functions A^. 



i2 roo 

Bi = ^ dqq^go{q)[Af\ql,R/l)f 
B2^^ dqq''go{q){Af\ql,R/l)f 



(5.9) 



These integrals are well behaved for finite R/l and the major contribution comes from the 
region in which the mesoscopic calculation is valid. This therefore resolves the difficulty 
encountered with the point scatterer, but the result is only valid for small R/l. 



6. Conclusion 

In the mesoscopic regime, we have shown how to calculate the effective charge and 
dipole moment that would be seen in the change of the transmitted light due to some 
simple objects placed inside a diffuse medium. The simplest calculation was for a single 
point scattering object. In this case we were able to sum up all higher order terms in 
the pertubative expansion (in the strength of the vertex) which contributed to leading 
order in the mesoscopic approximation. However, the point nature of the scatterer was 
incompatible with the mesoscopic approximation and this was reflected in the divergent 
result. To resolve this problem we considered extended objects. A collection of point 
scatterers gave finite results, but it was not possible to analytically sum all the higher 
order terms. We therefore turned to a spherical object and our main results are the 
formulae ( ^.8| ) for the charge and dipole moment of a small sphere in a multiple scattering 
medium. These results include an infinite sum of higher order corrections and despite the 
many simplifying assumptions that went into their derivation, they should be useful as an 
independent check of the numerical methods in common use. 
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Appendix A. Notation 

In this appendix we recall some basic definitions. The background scatterers that 
compose the diffuse medium have a density n and are individually described by the t-matrix 
t. The mean free path I is related to the elastic cross section by I = l/na^i = 47r/ntt. 
The green function for amplitude propagation in this medium has its wavevector shifted 
away from the value = cu^eo of the free scalar wave equation, V^V + ^^V' = 0) 
obtain Go(p) = {p^ — k'^ — nt)~^ = (p^ — k% — ikii/lex)~^ ■ In the second form we have 
explicitly split the t-matrix into its real and imaginary parts. The real part acts to shift 
the frequency, kj^ = k'^ + n^et. This effect is ordinarily described by an index of refraction, 
we shall take this as understood and drop the subscript on k. The small imaginary part 
of the t-matrix is most important and is related by the optical theorem to the total cross 
section including both the effect of elastic scattering and absorption, this is described by 
the extinction mean free path, n^mt/k = 1/lex- Absorption is traditionally described by 
an albedo factor that relates the absorption and elastic cross section by a^i = aatot or 
o'abs = (l/o — l)cre/ and also the mean free paths, lex = oZ. We shall only allow absorption 
by the extra, not background, scatterers. The optical theorem and relation to the various 
mean free paths and albedo reads, l/l = nti/ATr = an^mt/k = a/l^x- 

Appendix B. Integrals 

The following pair of integrals, related to the diagrams in figure 1, appear throughout 
the work: 



A mesoscopic approximation to these functions is sufficient for our purposes. 

The first integral /i is the ladder kernel of the SM equation and may be obtained 
most conveniently from the mesoscopic real space form as the Fourier transform of A^o = 




(B.l) 



e-^/747rZr2, 




(B.2) 



Hence the free green function of the SM equation is ^0(9) = 1/(1 ~ ^i(90)- 
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The integrals ([B.l| ) can also be evaluated directly by performing the integrals over the 
magnitude of the loop momenta, and then expanding the angular integrals in 1/kl to find 
the following mesoscopic limits: 

/l(qi)-5^-^l(«li) 

% (B.3) 

^2(qi,q2) ^A2(qi,q2) 

The properties of the angular integrals Ai,A2 are discussed in appendix C. 



Appendix C. Angular Integrals 

The definition and properties of the angular integrals Ai, A2 are given in this appendix. 

A / \ 1 /" 1 tan~^(q') ,^ ^ , 

M{q) = — dp — r = — C.l 

47r Js, 1 + ^q-P q 

The integral Ai is real and depends only on the length of the vector q. Note that Ai{0) = 1, 
and that the expansion for small q is Ai{q) = 1 — + g^/S + 0{q^). At large q, Ai 
behaves as 7r/2g. 

We also require a slightly more complicated integral, 

^2(qi,q2) = ^ [ dp J . ^ . (C.2) 
4:7T Js^ 1 + zqi.p l-zq2.p 

This expression is again real due to the symmetry p —>■ —p. It obeys the following, 

A2(q,0) =A2(0,q) =Ai(q) (C.3) 

Although a simple expression for this integral in not available, for the purposes here 
it is sufficient to consider an expansion in the incoming momentum. To leading order in 
small qi we have, 

A2(qi, qs) ~ A^{q2) - ^^(l - Aife)) = ^ife) - qi.q2A3(?2) (C.4) 

'?2 

where we have introduced As{q) = {l—Ai{q))/q'^. If both gi and q2 are small, A2(qi, q2) ~ 
1 — qi.q2/3. Some further properties involving angular integrals of A2 are also needed. 



/ dq A2{q,q2)=Ai{q)Ai{q2) 
, (C.5) 
dq qi.qA2(q,q2) = qi.ci2q'^ As{q)A3{q2 
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(D.l) 



Appendix D. Matrix Formalism 

Consider each of the 3 diagrams in figure 1 as a separate entry in a 3 x 3 matrix. The 
matrix version of the vertex for small incoming momentum (3.7) is 

Vi = -^QiAi(p'0 +P-P'47r/2PiA3(p'/) 
Qi = 7Z7 I ^V2A: 

-it/2k 



Pi = 7Z I ^V2/c 

-it/2k 

This form is preserved at higher order, though Q and P are no longer diagonal, and the 
measured charge/dipole moment is given by the sum of all entries of the respective matrix. 

The advantage of the matrix formalism is in connecting vertices using a propagator 
that can take account of the appropriate combinatorics discussed in the text. To only 
include ladder diagrams we write, 

G= 1 1 I (D.2) 

To also include crossed diagrams we have, 

G= 1^2 2 1 I (D.3) 
The higher order diagrams are summed to give. 




Q = Qi(l + — GQi 



(D.4) 



When evaluated using the ladder propagator we recover the result (3.12). When the crossed 
diagrams are included we must explicitly drop terms that are subleading in the mesoscopic 
approximation and are left with the result, 

Q = o — + Qi^(l/ae - l)Bi = ^ + —QjBi (D 5) 

+Pi3aez/(l/ae)S2 = — — +47rPfS2 (D.6) 



1 + 6(7ez/(l/ae)S2 ^ ^ ' - ^ ^ STrPi^a 
where the added terms correct for the fact that the two-scatterer diagram has no time- 
reversed counterpart. The main change is thus a factor of two in the S's compared with 
the ladder result ( |3T2D . 
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Appendix E. Classical scattering from a sphere 

We collect certain results concerning amplitude scattering from a sphere . A scatter- 
ing solution for a field obeying the wave equation, with incoming part e*^^, and symmetric 
object, is 

= Y,'WM'2l+l){ji{kr) + aihl'\kr)^ Yio{9) (E.l) 
I 

The coefficients ai are determined by the boundary conditions. Besides Dirichlet (D) and 
Neumann (A^) conditions, a more general possibility is to assume a surface impedance, Z, 
and require, dfcj) = ikcj) = Zdr(f>, on the surface. These coefficients are the same as the ones 



appearing in the green function (|5.1 



By calculating the fluxes we flnd the following cross sections. 



O'abs ~ 



A;2 

I 

Ctot = CTel + CTabs 

For D and conditions, aabs = 0, and for a large sphere {kR ^ 1), atot = Syrii^ 



Appendix F. Integrals for Spherical Object 

In the text we have only discussed the position space form of the bare vertex for 
the spherical object. The following representation in momentum space is convenient for 
calculations, 

yi(p, p') = J d^r J rfVe-^P-''e'P'-''Vi(r,r') 

'^'^ [ ^ 2 /" , / i2Sr^(- (O)-(l) , -(0) (1) , , - (1)-(1)\ 

= — drr drr y^K^^i'di 9v + aigy gl ' + +aiai'gl ' gy ) (p.l) 

^ex J R JR 

X [dr [ dr'e-'P-'e'P'-'-'^^^^-^^^^-^Pi{cos^)Pi,{cosj) 
J J A-K Atx 

Where 7 is the angle between the two vectors f and f'. We have written the mean free path 
as lex in this formula to prevent confusion with the angular momentum quantum number. 
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Although not tractable, this form can be used to prove the following two formulae which 
are the analogs of the integration formulae ( |C.5|) in appendix C. 



d^r j dVa(r)6(r')e^P'-''Vi(r,r') 



3 JR Jr P'r' 



j d^r j rfVa(r)6(r') z(p.r)e*P'-'"Vi(r,r) 



P P 



(F.2) 

Where a(r) and h{r) are smooth functions. The asymptotic form of the Bessel functions 
has been used and oscillating terms which contribute at higher order in in 1/ kl have been 
dropped in accordance with the mesoscopic approximation. 
The functions A^^^ and Ag^'* are defined as follows, 

A[^\ql,R/l)= r^^jR^e^R-r)/i = l_ f dp^^^ 
Jr I qr ^T^Jso i+«q-p 



4:TT 1+iq.p/ 
A^^'Ual B/l) - ^ r — ( ^ - e(«-^)/^ 



(F.3) 



They have the following limits as <^ /, 



A[''\ql,R/l)^A^{ql) 

^3 {.qlR/i) ^M{qi) = — j-^^ — 
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